
NPS ARCHIVE 
1964 

WILLIAMS, R. 




DU01.EY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOi 
MONTEREY CA 93943-5101 



Library 

U. S. Naval Postgraduate School 
Monterey, California 



This document has been approved for public 
release and sale; its distribution is unlimited* 



// 




DEMAND DURING LEAD TIME 
IN AN INVENTORY CONTROL SYSTEM 



a & # & a # 










# 



>4 w 



4F 

v*4* 



# 



# 

* 



* 



DEMAND DURING LEAD TIME 



IN AN INVENTORY CONTROL SYSTEM 



by 

Robert L. Williams 

0 

Lieutenant, Supply Corps, United States Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 

IN 

OPERATIONS RESEARCH 

United States Naval Postgraduate School 
Monterey, California 



19 6 4 



IU4 



(CX'/V\ S> , K , 




% 



Library 

U. 3. Nava] Portgraduate School 
Montc.ey* California 

DEMAND DURING LEAD TIME 

IN AN INVENTORY CONTROL SYSTEM 

by 

Robert Lo Williams 



This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

OPERATIONS RESEARCH 
from the 



United States Naval Postgraduate School 



ABSTRACT 



Demand during lead time is considered as the sum of a random number 
of random variables. Discrete random variables are considered as de- 
scribing demand and lead time and probability generating functions are 
used to determine the probability generating function of demand during 
lead time. The probability distribution of demand during lead time is 
extracted from the probability generating function by using a power 
series expansion and the calculation of probabilities from the expansion 
are discussed. The mean and variance of demand during lead time are 
also obtained from the probability generating function. 
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Introduction 



The two major problems in any inventory control system are to 
decide upon the amount to order each time replenishment is necessary 
and to decide what should be the minimum stock levels, of any items, 
at which an order should be placed. The second problem is the subject 
of this paper. 

The minimum stock level at which an order for replenishment of an 
item is the reorder point, denoted by R. The reorder point is a func- 
tion of the demand and the lead time of an item* The demand for an 
item is the rate at which the stock of the item is depleted. The lead 
time of an item is the length of the interval from the time a replenish- 
ment order is placed until the order is in stock, ready to satisfy the 
customer's demand. 

A simple deterministic model of an inventory system with constant 
demand and lead time is shown in figure 1. 
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Deterministic Inventory Model 
Figure 1 
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The stock is depleted at a uniform rate from a stock level of S 
until it reaches the reorder point. When the stock level of R units 
is reached, a replenishment order for S units is released and the 
order is received at time A, just as the stock is exhausted. The 
stock replenishment raises the stock level back to S units and the 
process is repeated. 

Rarely is demand uniform, however, and there is uncertainty about 
when the R units of stock on hand at the time of placing an order will 
be exhausted. A more realistic model of an inventory system with 
variable demand and constant lead time is shown in figure 2. 
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Variable Demand Inventory Model 



Figure 2 



S, R and A are the same in figure 2 as they were in figure 1. 
The expected depletion of stock is along the line S A. If the lead 
time was deterministic, an order placed at R, time X, would be re- 
ceived at time A. In the time interval during which replacement 
stock is an order, X A, there can be a variation In the demand j 
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therefore, the demand may deplete the R units of stock along any line 
between a maximum demand, R A*, and a minimum demand, R A”, shown by 
the shaded area A* R A". The variation in demand allows for the pos~ 
sibility of depleting the stock at time A* and not being able to meet 
demands for stock in the period A' A. There is also the possibility 
that the demand could be as low as R A" and in this case there would 
be B units of stock on hand when the replacement order is received * 

A third and the most realistic model, considering both variable 
demand and lead time, is shown in figure 3. The first two models 
were special cases of this model and therefore model 3 will serve to 
describe the inventory system considered throughout the rest of this 



paper. 
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Variable Demand and Lead Time Model 



Figure 3 



The variation in the lead time will cause the replenishment 
time interval X Y to vary. Time X is the time at which the reorder 
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point, R, is reached, but it is no longer certain that the order 
will be received at time Y„ The order may be received at time U or 
time V or any time in between, if we let time U and time V represent 
the minimum and maximum times of receipt of the order . Figure 3 is 
similar to figure 2, but the zero stock level of figure 2 is the 
same number of units below R as C is in figure 3. Figure 3 then is 
the same diagram except that the entire stock is Increased by C units 
to allow for the variation in lead time as well as the variation in 
demand. Let us consider the four possible extreme cases of maximum 
and minimum demand and lead time. A replenishment order is placed 
when the stock on hand is R units at time X. If the maximum demand 
and the minimum lead time occur, the stock will be depleted along 
line R A' until there are D units at time U when the replacement stock 
is received. If the minimum demand and the maximum lead time occur 9 
the stock will be depleted along line R A" until there are E units in 
stock at time V when the order is received. Minimum demand and mini- 
mum lead time will result in F units of stock on hand at time U when 
the order is received. In each of these three cases, R units of stock 
have been more than enough to provide protection against stockouts 
with both variable demand and lead time. Only in the case of maximum 
demand and minimum lead time would the stock be below C units. The 
last case is the combination of maximum demand and maximum lead time. 
The stock will be depleted along the maximum demand line R A 11 and will 
be exhausted at time Z. The order will not be received until the maxi- 
mum lead time, time V, therefore the system will be out of stock for & 
period of V-Z time units. 
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From figure 3? it can be seen that if the demand follows the 
line R A and time Y is the lead time* then the demand during lead 
time is equal to R - C units of stock. It can also be seen in 
figure 3 that as the height 0 C is increased the intersection of 
the actual demand and the actual lead time, whatever they may be, 
will intersect above the zero stock level more often than below. 

The reorder point, R, is then the sum of the expected demand during 
lead time plus a safety level or buffer stock of C units. The re- 
order point equation is thenj 

R = Z + V S L (1.1) 

where 

R = reorder point 

Z = expected demand during lead time 

V S L = variable safety level 

The safety level is denoted as a variable safety level because the 
protection against a stockout increases as the number C is increased. 
The safety level can be varied so that no stockouts occur or some 
small number of stockouts occur in a specified number of replenish- 
ment periods. That is, the inventory manager may decide that a 
safety stock large enough to protect against no stockouts would be 
too expensive and therefore he is willing only to maintain a safety 
stock that will protect against stockouts ninety or ninety-five per- 
cent of the time. The safety level is also a function of the reorder 
quantity since the larger the reorder quantity for a given lead time 
the less frequently the system is exposed to a stockout. This 
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consideration is taken into account when the protection level is 
determined and does not affect the reorder point equation J 9 > 1L 
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2. Discussion of the Problem 

The availability of alternative courses of action implies the 
need for managers to make decisions » The selection of the reorder 
point for an item in the manager's stock is then a decision problem 0 
The reorder point rule must satisfy an objective function (such as 
maximum service) subject to some constraints imposed by higher 
authority, available funds, available facilities, etc* In almost 
all cases, the method by which the problem is solved requires a model 
of the operations * Such a model is based on a set of relations link- 
ing the important variables,. Figure 3 is the geometric model of th® 
inventory system being considered and the reorder point equation (Id) 
will lead to the mathematical model which will be used to determine 
the probability distribution function of demand during lead time* 

The reorder point equation is a rule by which the inventory 
manager may decide when to reorder replacement stock* The variable 
Z is the mean of the probability distribution of the demand during 
lead time and the variable V S L is a function of the variance of the 
same probability distribution Therefore! in order to use the reorder 
point equation, the probability distribution function of demand during 
lead time or at least the mean and variance of the distribution func- 
tion must be known* 

Let Z be a random variable which describes the demand during 
lead time, then the number of demands that occur during a lead time 
may be written ass 



Z - D.j + I>2 + »«» + D^ (2*1) 
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where 



Z = the number of items demanded during a lead time 
D^— demands , i = 1, ••• > L 
L = lead time 

It is assumed that the demands are independent and identically dis- 
tributed and demand and lead time are independent discrete random 
variables with known probability mass functions „ In general* this 
paper will only consider discrete random variables since demands 
naturally occur at discrete intervals and lead times are normally 
described in days* weeks or months. Equation (2gl) is then a model 
which describes the demand during lead time by linking together the 
random variables D and L and is a model that can be used to deter- 
mine the expected demand during lead time* Z* and the variable 
safety level*, V S L, Equation (2d) can be described more concise- 
ly as s 

L 

y 

Z = »i (2-2) 

Using equation (2.2) as the model of demand during lead time* 
it will be shown that the mean and variance can be computed for any 
probability mass functions that describe demand and lead time or 
from the mean and variance of demand and lead time if that is all 
that is known about them. It will also be shown that the distri- 
bution of demand during lead time can be found* in the discrete 
case* by the use of probability generating functions. After the 
probability distribution is found* it will then be shown that it 
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can be used to set the reorder point using a variable safety level; 
given the level of customer service desired „ The final section 
treats a special case in which lead time is a continuous random 
variable and the resulting distribution of demand during lead time 
has a familiar negative binomial form 0 
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3* The Mean and Variance of Demand During Lead Time 

The mean and variance of Z are derived using general probability 
mass functions for D and L« 

Lets 

p (d) be the mass function of demand with mean D and variance 

p (1) be the mass function of lead time with mean L and vari- 
ance CT 

p (z) be the mass function of demand during lead time with 
mean Z and variance 

The mean of Z is the expected value of Z and is denoted as Z 

E(Z) = E [e(z)L)] 

= E [l E(D)] 

= E(D) E(L) 

= D L 

Therefore the expected demand during lead time is 

Z = D L 

The variance of demand during lead time is derived as follows s 
Var(Z) = E [var(z|L)] + Var f*E(Z|L)] 

Var(Z) = E [l + Var [l ECO)] 
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= E(L) (T d 2 + [E(D)j 2 Cr L 2 

- s - r - 2 _2 -.2 

= l cr + d o 

D L 

The mean and variance of demand during lead time are functions 
of the means and variances of the demand and lead time only, there- 
fore a knowledge of the mass functions of demand and lead time is 
sufficient to calculate the mean and variance of demand during lead 
time. They may also be calculated empirically using historical de- 
mand and lead time data to estimate the means and variances of demand 
and lead time. The mean and variance of a distribution function may 
also be calculated from the probability generating function of the 
distribution. The probability generating function of a probability 
distribution function is discussed in the next section and the calcu- 
lation of the mean and variance from a probability generating func- 
tion is shown in Appendix I. 
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4. Distribution of Demand During Lead Tima 



In this section, the probability generating function of demand 
during lead time will be derived in general* The generating func- 
tion will then be used to obtain the probability distribution of 
lead time demand by expanding the generating function in a power 
series and extracting the coefficients » Probability generating 
functions are defined for discrete random variables only* There- 
fore; the discussion will be limited to discrete random variables. 
The probability generating function of demand during lead time 
is the generating function of the sum of a random number of random 
variables and is defined as 



To derive h(s), let the probability generating function of demand 
and lead time be 




g(s) = Eta®*) 



f(») = E(s D ) 



and 



h(s) = E(s Z ) = E (s Z ]l = 1) 



now 
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since the Dj_ were assumed independent and identically distributed. 
Also since f(s) = E(s^), we have 



E(s Z |l = 1) = [f(s)] 1 

and 

E(s Z |L) = [f(a)J L 
E [e(s Z |l)] = E|[f(s)] L | 
= g [f(s)] 

since 

g(s) = E(s L ) 
h(s) = g [,(s>] 



Thus, if the distribution functions of demand and lead time are known 
or assumed, the probability generating functions of demand during lead 
time can be derived. Then with a knowledge of h(s), the probability 
distribution of demand during lead time is known. Examples showing 
the computation of the distribution of Z will be given later. 

Implementation of the model given in equation (2.2) may generally 
take place using familiar discrete distributions for the variables, 
since such distributions have convenient mathematical properties. 

These familiar distributions usually have positive probability at zero 
and this situation usually is not compatible with lead time, Prac- 
tically, lead time can never be zero, if for no other reason than it 
takes some time to prepare an order. On the other hand, if the data 
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in the inventory control system is based on the month as the basic 
time unit;, a lead time of one day could be considered zero. In any 
event, if the mean of the lead time distribution is small, then 'the 
use of a popular distribution will result in probability at zero which 
is not negligible. Mathematically, the model is correct since when L 
equals zero the sum is empty and Z is zero. Demand during lead time 
can also be zero when demand is zero during the lead time or when 
both demand and lead time are zero. 

There are several methods by which the input distribution can be 
adapted to consider either allowing or not allowing some probability 
for lead time to be zero. In the first case, when lead time can be 
zero, the straight forward use of the model with distribution ftme= 
tions of demand and lead time that are defined on the non-negative 
integers will have provided the desired results. There are two methods 
that may be used to prevent the consideration of zero lead time in the 
model. These methods redistribute the probability mass functions from 
a range over the non-negative integers to a range over the positive 
integers. The first is translation, i.e., let L = 1+X, where X takes 
non-negative integral values? and the second is truncation, 1 ©., 

use a frequency function that is conditioned on the variable being 
positive. Let us consider the Poisson distribution as an example. 

The Poisson distribution is usually defined ass 

-A x x 

p(x; A ) = - — , for A>0, x * 0,1,2,... 

xl 

The mass function can be translated over the positive integers by 
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writing it ass 



— ^ X*“1 

p(x| X) = ® — , for A > 0, X = 1,2, . 0 . 

(x-1 ) ! 



or the distribution function can be truncated by using a conditional 
probability,, That iss 



P fx = x|x > o] = 6 X A * - 
L J .1(1-.-*) 



\ X 

i A 

e A -1 x! 



A x 

o © p(x “ , f or A ^ 0, x 1, 2,ooo 

(e*-l) x! 



There are then three choices offered for adjusting the input to fit 
the situation or the preference of the user. 

Let us now consider some generating functions of demand and lead 
time so that they will be available for use in determining the prob- 
ability generating function of Z. The Poisson and negative binomial 
distributions have frequently served as good approximations for the 
distribution of demand so they will be used in the following examples. 
Lead time distributions have not been as well documented, but the 
negative binomial distribution is a good model for discrete waiting 
time random variables and when the parameter r equals one, the nega- 
tive binomial becomes the geometric distribution. As a second lead 
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time distribution, the Poisson distribution will be used* The use of 
the Poisson distribution for lead time may not be a good example , but 
with a large mean, the Poisson distribution roughly approximates a 
normal distribution with the mean equal to the variance fl^jo 

The generating functions for the Poisson and negative binomial 
distributions in each of the three cases are given below „ 

For the Poisson, when 



p(xj A ) = 



-A v x 

e A 



, X — 0,1, ooo 



The probability generating function iss 



CO 



f ( s) = E(s x ) = l \ s k p[x - kl 
k=0 L J 



- e 



X(s=-1 ) 



(4.1) 



when 



p(x| A ) = 



- A s x-1 
e A 



(x-1)! 



, X — 1 , 2 , ooo 



f(s) = s e 



A(s“1 ) 



(4.2) 



when 



p(x$ A ) = 



A 



( -1 ) x ! 



Xs „ 

f(s) = -4 - - 



e* -1 



, X * 1, 2, ooo 



for negative binomial distribution, when 
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p(xjr,p) = ^ r+x “^ P r q X 



» x — 0, 1 , 



f(s) 



= M-) r 

\1-q s J 



( 4 , 3 ) 



when the distribution is translated 



p(x;r,p) = ( r+ *^) f q*' 1 



> x — 1 , 2 s 



f(s) = 



■ . H-Y 

Vl-qs / 



when the distribution is truncated 



p(x$r,p) = 



hr 1 ) 



r x 

p q 



i-p 



* x = 1 , 2 , 



f(s) = 



_ p r [(l-qs)~ r ~l] 



1"P 



As an example of determining the probability generating func~ 
tion of demand during lead time, let 

D be a Poisson random variable 
L be a translated negative binomial 

then 

f ( s) = e A(s “ l) 

and 

h(s) = g jjf (s)J 
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r — i 






= p r 


1 e*< 


s-l)k 


/r+k=2\ q k--1 




k=1 




\ k-V 






(DO 

SP 




U A 

II 


e A(s-1 ) 


X 

k=1 


e A(s-1) (k-1) ^r+k~2j 






00 

V"~l 




II 


e A(s-1) 


2 


/r+k-s'j [ qe AU-l)]*" 1 






k-1 




= F r 


e A(s=l) 


jj-q© 


A < s— 1 ) ] -r 



The probability generating functions of demand during lead time,, Z 
were derived since they can be used to find the probability of d®=> 
mand during lead time., That is, sine® 

co 



h(s) = E(s Z ) 



s k P 



k=0 



9 



Then in the power series expansion of h(s), the coefficient of 
is the 'M- For example, writing out several terms of h(s), 
we haves 



+ s P £ * 1 

It therefore is necessary to put h(s) in the form of a power series 
in order to determine the probability of demand during lead time 
for any k. 

Using the probability generating functions for the variable Z 




] + s 2 P jz - 2 
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in the previous examples, we can put them in a more useful form* 
For Poisson demand and translated negative binomial lead time 



h(s) = e A(s ' 






= e A(s “ l) p r 



I hi 



= P r 



= P r 



oo 

z ( r+ !‘- 1 ') 

k=0 V k / 



q k 3 s PS (k+l) - (k+1 ) 



OO OO 

Y ( r+k-l\ q k e - A (k+1) Y s n 

k=0 \ k / n=0 n! 



CO oo 



= Z s “ P r z ( r+k - i 'i 

n=0 k=0 V k ; 



q k e - A(k+1) 



n! 



and 



© 

o e 




is the coefficient of s n in the power series 



oo 



& - ■] ■ > r £ hi 



q k e - A(k+1 ) CA(k+l)] n 

n! 



co 



- rJ 

- p e 



£Z. 

k=0 V 



n! 



^r+k- 
k 



’) (qe A ) k (k«) n 



The calculation of these probabilities can be reduced if the sum 
from zero to infinity is reduced to a finite number of terms with a 
very small remainder, A method for finding the number of terms to 
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be summed until the remainder is negligible will be shown along with 
a computing technique for determining the probabilities. 



To find a bound for the remainder of the series, let 




(qe-* ) k (k+l) n 



a - qe 



and 



co 



Ri 



y y u 

^ u . = %+1 ^ — 



K+1+1 



k^+1 Uk UR+1 1=0 U K+1 



then 



U K-H+1 - (k+1+l) a k+1+1 ( k+2+l) n 

Vi (^) a K+1 ( k+2 ) n 



(K+1+1) !(r+K) 



/ K+2+l' 

V K+2 






= (r+K+1) (r+K+1-1),., (K+l+r-(r-2) ) 

(K+r) (K+r-1)... (K+r-(r-2) ) 



/ K+2+l \ n 
\ K+2 / 



a 



1 



(r+k+l) (r+K+1-1 ) . . . (K+l+2) / K+l+2 ^ n & 1 

(K+r)(K+r~1 ) ... (K+2) V K+2 J 



Dividing term by term 
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then co 

r k ^ “k+i ^ ^ + *y) &1 

It is possible t© dominate the right hand side with a integral pro- 
vided the integrand is monotone decreasing. This will be true if K 
is sufficiently large. Thus, let 

-B s in a = = A + In q 

and 

¥ (1) - (i + 4=) n+1 .- B1 

then 

T ' (1)= ( 1+ t )" # ' b1 [-( 1 + t) b + t ] 



Setting the derivative equal to zero, we have 




B 



o n+1 

. * . the maximum for 1 occurs if K > — - 

B 

hence we must us® K such that it is at least this large. 

New f or K > 

B 
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CD 





+ 



+ 



-i.) n+1 .1 

K / 

1 \ n+1 -Bx 

TV e 



dx 



= 1 + 





dx 



making a change of variables 9 let v=Bx, then dv=Bdx, and 



n+1 ~4E> 

2 o + +) al 1 + 2 ft 1 ) 0 I vi e " T 



dv 



n+1 



= 1 + 



2 ft 1 ) 0 1 s 



& i + 



(n+1) I 
B 



(, + -If 1 

V KB/ 



which is a bounded function of K„ thus 



„ r (n+ 1 )! I 1\n+n 

s k«“k+i L ,+ — rs) J 



a — d-t* {■*=)"] 



now if KB > n+1 
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%4(w) *“ B(K+1> (*«)' [i + .] 



In practice^ all of the constants and parameters in the above 
equation are known except K„ K can be determined by setting the 
equation equal to a very small number and solving for K, The 

error in P (Z~n) will then be less than that small number times 

y =1 1 n 

P e . 

n! 

New that the value of K has been approximated;, the probability 
distribution of demand during lead time can be found as 



K 



[i . sgi («>-■) 



a k (k+1) n 



In order to simplify the calculation of this number manually or in 
a computer s the technique of "nesting** will be shown,, 

K 






k=0 



r+k - 1 ' a k (k+1) n = 1+(^ aa^^ 1 ) 



2 -n , 

a 3 + o o « + 



( r+k - 1 )a K <K + ,)“ 



1 + ra 2 n + a 2 3" + ... + a K (K+1) n 

21 ft I 



- |l + ra 2 n ^1+ 



(r+1) 3 D 



))•")) 



This can then be written in the compact form 0 
Let 

r+K=1 / K+1 \ n 



K 



a 



(k* ) 
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and 




B i = 1 + Ai B q 



o 




where 



n 




and B ^ = 1 



then 



P [z = n] = (j 




As another example 

Example 2 

Let 

D be a negative binomial random variable 
L be a translated Poisson random variable 



Then 




from eq 9 (4.3) 



and 




from eq. (4.2) 



Then 



h(s) = f (s) g [f(s)] = f (s ) e , X C f(s) - 1 3 
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Let us again find the Maclaurin expansion 



h(s) = 



tester -■]{ 



CO 



_ a A r ( k+1 ) / i \ ~r ( k+ 1 ) 

- e “T? p U-qs; 

k=0 k! 



oo co 

A ^ A r(kH) 2 /r(k+l)+n-l\ / \n 

k=0 “Id P n=0 V n J (q s) 



<x> 



00 



1 » n e- A q n 2 ( r(k+1)+n - 1 ) p r (k+l) 
nsO k=0 ' n ' k! 



and 



CD 



i* <■ p r g ( r(k+ : )+ i 



(Ap r ) k 
k! 



To find a bound for the remainder of the series 



Let 



K K 



*K 



= I n k = 1 /r(K+1 )+n-l\ Oi£? 

k=0 k=0 V n / k 



and 



% 



£50 O t> 

y . y u k+i 

■ . 4* . U k ^ U K 



k=K+1 



Then 



_ fr(k+l)+n~l] ! fl k 
' k (r(k+l)-l] ! n! 



, where a = Ap r 
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and 



U K+1 (r K+rl + n-l)! (rK + r-1 ) 1 i 
- a 



u; 



K (p K + pl-1 ) ! (p K + r+n-1 ) 1 



(rK+ r!+r+n°1 ) (rK+rl+n-2) , «, «, (rK+r+n) & 1 
( rK+r 1+r “1 ) ( rK+r 1+ r-2 ) . . , ( rK+r ) 






1 + 



n ) 


L , n ) ( 


'l+ ° ) 


p(K+l+l)-1 ' 


I 1 ' J® © 0 l 

' r(K+l+1 )-2' ' 


r ( K+1 ) J 



Thus 



< (’ * -ass/ •' 



U, 



K+l 



r -i 1 



.,y •] 



let 






We then have to choose minimum K such that r<^ so the resulting 
geometric series can b© summed, i.e„, 



*• [(■ * afef •] < ’ 
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if 



1 + 



n s' 

Hk+iT < 




and 



K> 




1 



Then IS ■*£ 1 
Now 

00 <x> 

o y u k+i y y y i 

K = u v /L, -T- < U 0 

K K 1=0 K K i=o 

1 

= u v — — 

K 1~ Y 

and to solve for K we have 

-X n r ( 1 ^ _ 0 (K+l)+n-l] 1 p rK e“ A q n p r 

q P UK ll = V J [ r (k+1) -l3 ! n ! [l- p* (1+ Jjj ) r ] 

Upon setting this equation equal to an allowable error, K can be 
determined o 

The estimate of the remainder is the value of the allowable 
error from which K was determined. 

Then, given K 

r(k+l)+n-l\ (Xp r ) k 

n I k! 
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The ^nesting'* technique is helpful in computing the probabil- 



ities of demand during lead time. 



P 




= e 



P r q n a k 



, r ^ 1 



and 



A - .( r+p ~'Q. .l + (2r+n-1 ) 1 (Ap r ) + (3r+n-1 ) ! (Ap r ) 2 + 
K (r-l)lnl (2r-l) !n! (3r-l)l n! 



= If+SzllL ( i + (2rfn°1 ) . . <..(r+n)_ ( A p r ) / 1 + . c 0 

(r“l)!n! V (2r-l) . (r) V 



o o + 



(rk+n-1) ... (rk+n-r) (Ap r ) 



(rk“1 ) 



(rk-r)k 



^ + O O 0 



1 + 



(rK+n-lJ ..o (rK +n-r) (A p r ) 



(rK-1) 



(rK-r) 



))•")) 



Thus, if 



B. 

J 






D* 



n 



r(K-j). 






and 



Cq = 1 , then 

c i = 1 + Bq c 0 



>1 " 1 + B j °j 
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and 



4k = (r-l)ln! (’ + ®K-2 °K-2 ) 

Now that the distribution of demand during lead time and methods 
for calculating the probabilities have been discussed, it is a simple 
matter to find the number of items of stock that are required as the 
reorder point. The reorder point must be the number of units of stock 
that will last until new replacement stocks arrive with the probabil- 
ity of a stockout during the lead time small. 

Let a = Acceptable probability of a stockout during lead time 



Then 

P [z>r] < a 

must be solved for the value of R. 

1 -P [z > r] = p [z = o] + F [z = l] + ... + P [z = r] 

.*. R is the smallest number such that 
R 
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5. Negative Binomial Distribution of Demand During Lead Time 

It would be very helpful if the distribution function of de- 
mand during lead time could be found for any distribution functions 
of demand and lead time in a familiar form, but this is not possible. 
It can be done in some special cases. One known case is given belcw 
as an example. Consider deriving the mass function of demand during 
lead time where s 

(1) The demand is random and the number of demands d 
during a fixed interval of time t has a Poisson 
distribution. 

A a t / .\d 

p(d|at) = yd = 0,1,2,... 

d I 

where a is the average number of demands per unit time. 

(2) The lead time, t, has a gamma type probability density 
function 

-ut / .\k-1 

f (t;u,k) = | U y t>0 

k k 

with mean — - and variance 

u u 

Then the probability of exactly z demands during the lead time is 

-ut . .k-1 

ML, di 



-<X> 



p(*J 



k,Uy 



a) = 



-at 



(at)' 



z! 
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-CD 



a z u k 



z+k-1 -(a+u)t _ 

t e dt 



; ! ( k“1 ) 1 o 



Letting x = (a+u) t, then 



z k 
a u 



%CD 



p(zsk s u,a) = 7- — ; »xg j ; — — - 

y 1 * * (a+u) zTR zl(k-l)! 



x z+k-1 e -x dx 



z - k (k+z-1 ) ! 



a u 



(a+u) 24 ^ z!(k“l)! 



f k+z-l\ / u \ k / a \ s 

' Z ' \a+uJ \a+u/ 



the last step being valid if k ii an integer,, Thus the demand during 
lead time has a negative binomial distribution with 



mean = 



ak 



u 



and 



variance 



= (1 + 2) 

u V u / 
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6. Summary and Conclusions 

The distribution of demand during lead time has been discussed 
when discrete random variables describe the demand and lead time. 

It is unfortunate that the distributions do not appear as well 
known distributions. The use of probability generating functions 
has been shown as a method by which the distribution can be found 
and used. Calculations of the probabilities are cumbersome since 
in the examples shown each of coefficients of the in the prob- 
ability generating function involved an infinite sum. The amount 
of computer time required for the calculations has not been investi- 
gated, but since the calculations are straight forward and the number 
of terms required before the infinite sum converges to an approximate 
value are relatively small, the use of the procedure outlined in 
Section 4 should be useful and not too cumbersome to preclude their 
use. 

The use of the generating function method requires data for the 
estimation of the parameters in the demand and lead time distribu- 
tions. This data is readily available in the historical files of 
an inventory system. Other methods of estimating the demand during 
lead time requires collection of data directly from the number of 
demands during a lead time. This data is not usually collected and 
even if it were available sufficient data to obtain a good statis- 
tical sample would require a very long time. 

Further investigation of distributions calculated by the methods 
in Section 4 may show that the probability distributions could be fit 
to some known distributions and the resulting distribution used to 
approximate the demand during lead time. 
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This investigation has been made primarily to consider a method 
of determining the demand during lead time distributions for low 
demand items,, Low demand items are particularly troublesome in an 
inventory system. Military spare parts inventories have large 
numbers of very low demand items. The low demand items are also 
characterized by high variances £6], 

The use of the Poisson and negative binomial distributions 
have been used extensively in this paper since these distributions 
have been found to fit low and medium demand items. The Poisson 
distribution is generally used for items with a demand of two or 
less items per year. The negative binomial distribution has been 
used for items with demands ranging from two to one hundred per 
year [7], In the case where the Poisson demand distribution is 
applicable, the generating function method for establishing a re- 
order point can be readily applied. Since the reorder point will be 
relatively low, the number of probabilities that must be calculated 
before the sum is equal to or greater than the service level desired 
will be small. It has not been discussed in this paper, but when 
reorder points have been previously established, the probability 
generating function for demand during lead time can be used to find 
the tail of the distribution. With a reorder point previously estab- 
lished, a good estimate of the value of reorder point is available 
and with several trials on either side of the old reorder point, 
the tail of the distribution could be used to determine the new order 
point. In this manner, the probabilities would not have to be summed 
until the sum was greater than the desired probability. Feller [ 2 J 
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gives the following relations between the probability generating 
function and the generating function of the tail of the distri- 
bution. 

Let Q (s) = the generating function of the tail 



Then 

8 (s) = 

1 -s 

The calculation of the tail in this paper has not been investigated 
since it gets more complex than the probability generating function. 
However, in special cases and in a particular situation as described 
above, the tail of the distribution may be the easiest method for 
setting a reorder point. 

The special case of Poisson type demand and gamma distributed 
lead times was discussed because this combination of demand and lead 
time distributions lead to a negative binomial distribution of de- 
mand during lead time. Further investigation of this model with 
actual data may be very fruitful. The Poisson distribution is a 
good distribution for low demand items and the gamma distribution 
is a very versatile two parameter distribution that can be used to 
approximate a wide range of lead time patterns [5]]. Consequently, 
any demand and lead time process that may possibly be fit to this 
model should be considered, especially since the probability law of 
demand during lead time may be readily computed. 

The writer wishes to thank Professor Robert R. Read of the 
U.S. Naval Postgraduate School for his suggestions, assistance, en- 
couragement and cooperation in the preparation of this paper. 
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APPENDIX I 



MEAN AND VARIANCE FROM PROBABILITY GENERATING FUNCTIONS 
The mean and- variance of a discrete random variable may be 
calculated from the probability generating function. It may be 
shown that (see e.g., Feller [2]) 

E(X) = f*(l) 

Var(X) = f"(l) + f'(l) - [f • (1 )] 2 

As an example a the mean and variance of demand during lead 
time when lead time is a translated Poisson random variable 

A X_1 

p(xjA) - e“^ — - — , x = 1,2, 000 
(x— 1 ) ! 

with mean = (A +l) and variance = A 
and demand has the negative binomial distribution 

p(x$r,p) = ^ +X ^ p r q x , xs 0,1 ,2, 

with mean = and variance = ^ 

P . P 



then the random variable Z has a probability generating function 




taking the first derivative with respect to s we have 



and 



h'(l) 



rq 

P 



+ 



rq A 

P 




(l+A) 



the second derivative evaluated at one is 



h"(l) 



t rq(l + A )] 2 



2 2 
r q 



A + rq (l+A ) 



3 



then 



r- i2 

Var(Z) = h"(l) + h'(l) - [h'(l)J 

= (l+A) 2 + A + ^2 (1+A) + ^ (l+A) "^^j 2 ( 1+ A) 2 



These results can be verified against the general formulas for 
the mean and variance given in Section 2. 
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